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Section  1.  Introduction  ^ 

*r  > ' 

v  C  * 

;  For  independent  chi-square  variables^*,  and  x^/with  m  and  n  degrees  of  freedom, 
respectively,  we  consider  the  quadratic  form 

Q  =  «lXm  +C2xl 

where  the  positive  are  distinct. 

This  paper  gives  exact  finite  expressions  for  the  distribution  of  Q  in  terms  of  available 
functions  such  as  the  distribution  function  of  chi-square  random  variables,  modified  Bessel 
Functions,  Dawson’s  integral,  (tabled  in  Abramowitz  and  Stegun  (1964))  as  well  as  the 
distribution  of  Cj!)Xi  +  <42)Xi  (tabled  in  Solomon  (1960)).>These  formulas  are  useful  for 
checking  the  accuracy  of  approximations  and  tables  of  the  distribution  of  Q  and  provide  a 
simple  alternative  in  their  absence. 

For  large  m  and  n,  reasonable  approximations  to  the  distribution  of  Q  are  available. 
For  the  general  quadratic  form  Williams  (1984}^5ompares  algorithms  for  truncations  of 
infinite  series  expansions  of  the  distribution's  ee  Johnson  and  Kotz  (1970).)  Oman  and 
Zacks  (1981)  give  a  mixture  approximation  and  Davies  (1980)  provides  an  algorithm  for 
an  approximation.  For  small  values  of  m  and  n,  tables  for  the  distribution  of  Q  are  given 
by  Harter  (i960),  Johnson  and  Kotz  (1967),  Marsaglia  (1960),  Owen  (1962),  and  Solomon 
(1960). 

Distributions  of  the  form  Q  arise  in  a  number  of  applications.  Solomon  (1961)  noted 
that  probabilities  of  hitting  targets  frequently  reduce  to  the  distribution  of  quadratic  forms 
of  the  type  Q.  Pillai  and  Young  (1973)  show  that  the  trace  of  a  2-dimensional  Wishart 
matrix  is  distributed  as  Q  with  m  and  n  equal.  The  variable  arises  in  the  engineering 
literature  described  as  a  weighted  unbiased  Rayleigh  variate  of  dimension  two.  (See  Miller 
(1975)).  A  very  important  application  is  the  distribution  of  chi-square  goodness-of-fit  tests 
with  estimated  parameters.  Certain  two-sample  chi-square  tests  described  by  Moore  and 
Spruill  (1975)  have  asymptotic  distributions  of  the  form  Q.  Alvo,  Cabillo  and  Feigen 
(1982)  show  this  for  the  average  Kendall  tau  statistic.  The  distribution  of  Q  for  small  m 
and  n  for  the  average  Kendall  tau  statistic  is  provided  as  an  example  in  Section  3. 
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The  exact  expression*  for  the  distribution  of  Q  may  also  be  useful  for  approximations 
for  more  general  quadratic  forms,  especially  in  the  case  where  there  are  essentially  two 
groups  of  coefficients  nearly  alike  within  groups,  i.e.  the  distribution  of  Q  is  an  approxi¬ 
mation  for  the  distribution  of 

m-f  fi 

<3'=E 

«=1 

where  ss  clf *  =  1, . . . ,  m  and  a*  *s  ej, i  =  m  +  1, . . . ,  m  +  n.  The  exact  expressions  for 
the  distribution  function  of  Q  are  given  in  the  next  section. 

Section  2.  Exact  expressions  for  the  distribution  function  of  a  linear  combination  of 
two  chi-square  random  variables 

The  results  in  this  section  give  exact  expressions  for  the  distribution  function  of 

Q  =  «lXm  +  C2Xn 

where  the  positive  are  distinct  and  Xm  Xn  *re  independent  chi-square  random 
variables  with  m  and  n  degrees  of  freedom  respectively.  The  first  theorem  handles  the 
case  where  at  least  one  of  m  and  n  are  even.  Corollary  2.5  gives  an  expression  for  the 
distribution  of  Q  in  terms  of  that  of  a  quadratic  form  with  fewer  degrees  of  freedom.  This 
corollary  can  be  applied  repeatedly  to  give  the  distribution  function  of 

ciXjk+i  +  eaXa/+i 

in  terms  of  modified  Bessel  functions  Io  and  I\  and  the  distribution  function  of 

Vi  -  ci  Xi  +  Xi 

Tables  of  the  distribution  function  of  Qi  are  given  by  Solomon  (1960)  and  tables  of  I0  and 
1 1  are  given  in  Abramowitz  and  Stegun  (1964).  In  an  example,  a  representation  for  the 
distribution  function  of  e^xl  +  e^xl  “  given. 

The  following  theorem  gives  the  distribution  of  Q  unless  both  m  and  n  are  odd. 


Theorem  2. 


Let  Xm  *nd  xlk  independent  chi-square  variables  with  m  and  2k  degrees  of  freedom, 
respectively.  Then 


A 


p  VW  +  >  1  =  p  t*»  >  +  £ 0iP  [*«»-/>  >  «<,]  ■  n 


where 


0-(  Vf  £i -Y  Lilli 

1  \oo-oj  \|fll  —  Oo|/  r(f)(j 

If  ai  >  ao,~ii  is  P  [xm+jy  <  ai  ~  flo]>  If  oi  <  flo.  and  m  is  odd, 


**  ~  h  ^T) 


where  D(y)  is  Dawson’s  integral  tabled  in  Abramowitz  and  Stegun  (1964). 

Remark:  Note  that  the  result  in  the  theorem  is  completely  general  since  we  may  write 

P[cxxlk  +  c2Xl>c]=  + 

L  ®o  °i 

where  oq  =  cc^ 1  and  a!  =  ccj  1 . 


Proof: 


Hence, 


Equation  3.393,  #1,  p.  318,  of  Gradshteyn  and  Ryzhik  implies  that  the  integral  above  is 
0(y,j  +  l)iFi  (y.J  +  l  +  y;(ao-ai)/2)  where  XF\  is  the  confluent  hypergeometric 
function. 

For  ao  <  ai,  a  theorem  of  Bock,  Judge  and  Yancey  (1984)  implies  that  for  odd  m, 

r(y)  K°1  -  ao)/2]*  +  1  +  y;(ao  ~  <*i)/2) 

(-!)("-»)/»  <»(”x'>/»r(<-  j)(^7)*-|f  (xS< .,<«■-« o] 

(r(?  +  j  + 1))-'  (<  -  - 1)\ 

Applying  these  to  the  integral  we  have  we  may  write  (*)  as 

9  k-l 

e  ^at _ V^f— V 

r(?)(-2)=f1(«, 2 


r«  -  ^)(^)<-P  [xj,_,  <  -  ao] 

t  («  -  JSTil)!(j  +  -  0! 

Interchanging  the  orders  of  summation  and  setting  t  =  l—  above  gives  (*)  as 


lL«h»  v  E^ilkkl^zi  v  ilLP  (y» 
(5i  h  fcv-v  1 


X2i+m  <  «1  ~  «o] 


Because 


[*«*-•>  >4 


we  can  substitute  this  in  the  last  expression  for  (*)  and  the  theorem  is  shown  for  ao  <  a1 
and  odd  m.  A  corresponding  evaluation  of  jFi(y  ,j  +  l  +  y;  (do  -  «i)/2)  for  even  m  gives 
the  same  result  here  and  the  definitions  of  7,  and  complete  the  proof  of  the  result  for 
a,  >  oo  If  dj  <  do,  then  a  theorem  of  Bock,  Judge  and  Yancey  (1984)  implies  that  for 


S 
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Qn  =  W  +  coXn  where  Co  >  0.  For  c  >  0,  and  n  >  2, 


If  n  =  2, 


P  [W  +  Ml  >c]=  (2c0 )f9Jc)  +  P  [W  +  coxLa  >  *] 


P[W  +  cox?  >  c]  =  (2c0 )/Qs(e)  +  P  [W  >  c] . 


Corollary  2.3.  For  the  quadratic  form 

Q  =  ClXm  +C2Xn» 

we  have 

P\Q>e}  =  2 c2fQ{e)  +  P  +  c2xl-i  >  c) 

where  /q( x)  is  the  density  of  Q  and  Xo  =  0. 

Proof  of  Theorem  2.2.  Let  Qn  =  W  +  cox*-  Le*  ft J»(z)  *be  density  of  Qr 


/Q-(c)  =  ^[P[lV  +  tox’<c]] 


We  may  write 


r  „  ,  re,‘°  U  /-C-CoU  "1 

ptlv  +  -rcfP=75-  [/„  H* 


e  e-J/2(c-«)/co  r  rt 


=  risr 

/o  Co 


[/: 


r(f)2-/a  Uo 

where  t  =  c  -  cou  is  the  change  of  variable. 

Differentiating  this  last  expression  implies 

/■e  (£z!)$-3e-(c-t)/3e0 

*•  w  -  (2e°r‘</(" 2  3)  Jo  iy. 


dFjv 


dt 


/•e  /*e  (£^)n/J-le-(e-«)/aco  ,t 

+  /(ft  =  2)[/o  d*v]  +  (-i)/o  —  Cor(|j2va - l/o 


.  Then 


=  (2c0)_1{/(n  >  3)/>[W  +  coXn-a  <  c]  +  I(n  =  2)P[W  <  c]  -  P\W  +  c0xjl  <  «)} 


rt vyvvv 


;  ovrorKwvy 


The  following  theorem  gives  the  density  of  Q  in  terms  of  a  confluent  hypergeometric 
function. 

Theorem  2.4.  Let  m  and  n  be  positive  integers  and  let  ci,c2  and  c  be  positive.  Then  the 
density  of 

Q  =  ClXm  +  e2x5t 

is 

,  ,  ,  y(m+n)/a-le-*/ae.  ^  /fl  m  +  n  #  ^  _uy  ^ 

fQW)  r(=^)(2c1)'"/2(2c2)'‘/2l'Fl(2’  2  :(Cl  Cj  2 

for  y  >  0  where  \m  and  Xn  3X6  independent  chi-square  random  variables  and  1F1  is  the 
confluent  hypergeometric  function. 


Proof.  Let  W\  and  W2  be  independent  random  variables  such  that  Wy/cx  has  a  chi-square 
(m)  distribution  and  W2/e2  has  a  chi-square  (n)  distribution.  Then  the  density  of  Wy  is 


hx(x)  = 


zm/2— — */2c, 

r(f)(2Cl)-/2 


for  x  >  0. 


The  density  of  W2  is 


M*)  = 


zn/2-Ie-*/2c, 


for  i  >  0. 


Then  the  density  of  Q  =  Wx  +  W2  is 


/q( y)  =  [  Mv  -  *)ht{x)di 
Jo 


g-v/ici  (y  _  zJm/2-  lzn/2— lg  */2(e,  1  -e,  1  )J 

=  (2cl)"»/3(2c2)"/2r(f)r(f) 


The  integral  in  parentheses  can  be  written  as 


mnW"*"1’-'  F  l”  "  +  n  ,  ,  y 

r(  m+n)  1  l'  2  ’  2  ’  1  5  2 


c-»/2c,  (m+n)/2-l  f  I  n  m±n.  t  -l  _  -ls*\ 

fQ(y)  =  - - 1 - —LLiillZEjili _ >*>  II 

JQW)  r(=^)(2el)"*/a(2ea)'*/a  ,|1 

The  following  is  a  direct  result  of  Corollary  2.3  and  Theorem  2.4. 

Corollary  2.5.  Let  co,ct  and  c  be  positive  and  assume  that  chi-square  variables  are  inde¬ 
pendent  in  the  following  expressions.  Then  if  m  >  2, 

P\coXm  +  ClX,  >  C]  =  P  [coXm-3  +  clXl  >  c]  + 


For  m  =  2, 


c0(f)m+n/J-1<-e/JCo  p  ,n  m  +  n  c, 


*o  2  2  2 


PfcoXj  +  clXn  >  c]  =  P  [x*  >  ~]  + 

Cl  J 


£  +  1;  C-{c~ 1  -  cr1)). 


Remark:  Repeated  applications  of  Corollary  2.5  enable  one  to  evaluate  the  distribution  of 
Q  when  m  and  n  are  odd  since 

where  i*  the  modified  Bessel  function.  (See  Equation  13.6.3  of  Abramowitz  and 

Stegun  (1964).) 

Examples: 


(a)  For  c2  <  C\  and  y  =  |(c2  1  -  e1  *),  we  have 

P[eiXi  +  e3xl  >  cj  =  P  x?  >  — 

ci 

+  e^T(Jf-)ip(yi) 

»eiy 


P[ciX?  +  c2x<  >  c)  =  -P[x?  >  7*1+ 


Cl 


«*T(— )‘  '»(»*)(»  +  [1  +  V-1])  -  — ~-r 

*ciy  [  4e2  4c2y* 

(b)  For  d,  =  c/4ci,i  =  1,2, 

P[cix!  +  «2X?  >  c]  =  p[c\ i1}x?  +  c$3)x?  >  c]+ 


v/ 4di<i2  e  (dl+<i3^  {/o(d2  -  ^i)  +  fi(^2  _  d\)}  ■ 

^[cil)xl  +  ci2)xl  >-c]  =  P[c(i1}X?  +  4J)X?  >  «]+ 

e-(*+-.)(4rfl<i,)*  |/0(i1  -  d2)  +  -  <*a)  J  • 

For  instance  with  Cj  =  .25,  c2  =  .75  and  c  =  1.8,  we  get  after  substitution, 


t 


F(.25x^  +  .75*?  >  1.8]  =  .292. 


Furthermore  with  ct  =  ±,  c2  =  §  and  c  =  8,  we  get  after  substitution 


P 


1(1)  3  2(3>  - 
3  X3  +  3  X?>8 


.018318. 


It  is  instructive  to  test  an  approximation  for  the  two  probabilities  just  evaluated.  We 
replace  the  random  variable  Q  =  cix^  -t-  c2x„  by  [c\2p)k  and  obtain  values  of  c.p.  and  k 
by  equating  the  first  three  moments  around  the  origin  of  the  two  random  variables.  When 
the  three  parameters  are  computed  we  need  refer  only  to  the  usual  \2  values  for  p  degrees 
of  freedom  to  obtain  percentiles  of  the  distribution.  If  the  first  three  moments  of  Q  about 
the  origin  are  and  p'3,  we  have 


H  =  (2c)kT(k  +  v)/C 
p\  =  (2c)7kT(2k  +  v)/C 
p'3  =  (2c)3kT(2k  +  v)/C 


where  t;  =  p/2  and  C  =  r(r). 

It  is  convenient  to  define 

R2  =p'/p2  =  cT(2k  ^v)/{T(k  +  V)}2 
i?3  =  p'/p3  =cJr(3Hr)/{m-  +  r)}3 

We  first  calculate  R2  and  P3  from  the  moments  of  Q  and  solve  for  k  and  v.  then 
c  is  obtained.  Computer  routines  are  available  to  perform  these  operations  and  then  to 
calculate  probabilites  of  \2  even  with  non-integer  degrees  of  freedom.  Thus  probabilities 
for  Q  may  be  approximated,  since  we  have 

P{Q<<?}  =  P{X’  <(<J1/k)/c 

See  Solomon  and  Stephens  (1977.  19S0)  for  further  description  of  the  procedure. 

In  our  examples,  suppose  firstly  Q  =  .25\j  +  .75\j.  Fitting  {c\2p)k .  we  get  c  = 
0.13015, p  =  9.08011.  k  =  1.70.  and  so  Q  ~  (0.13015\$  080n  )]  70  and  P{Q  >  1.8)  =  .2925. 
By  exact  methods  we  obtained  P{Q  >  1.8}  =  .2920. 

Secondly  consider  Q  =  5X3  +  5X3-  By  equating  moments  we  get  c  =  0.32S51.p  = 
7.54926./;  =  1.18  and  Q  ~  (0.3285\~  1S  and  P{Q  >  8}  =  .01834  as  contrasted 

with  .018318  from  direct  calculations.  It  can  be  seen  that  the  ( c\7p)k  approximation  gives 
excellent  results.  Other  examples  are  given  in  Solomon  and  Stephens  (1977,  1980). 


Section  3.  Example:  the  average  Kendall  tau  statistic 

For  the  rankings  of  r  objects  by  n  judges,  the  average  Kendall  tau  statistic,  fn,  is 
the  average  of  Kendall’s  rank  correlation  between  each  of  the  (!J)  pairs  of  judges.  The 
null  hypothesis  is  that  the  r  rankings  of  the  judges  are  picked  at  random  from  a  uniform 
distribution  on  the  r!  possible  rankings.  As  n  — ►  00,  the  null  distribution  of 


(nf„  +  1)  • 


3r(r  —  1) 
2 


is  that  of 


11 


^*1 


V  N 

r.. 


^2*1; 


where  (r  j l)  is  the  binomial  coefficient.  (See  Alvo,  Cabilio  and  Feigein  (1982)  for  this  result 
and  discussion.)  The  results  in  this  section  are  derived  from  the  results  of  Section  2  using 
algebra. 

For  r  =  3, 

For  r  =  4, 

p  [s(1>xi  +  <’>x!  >  i]  =  P  [s<’>x?  +  <’>xi  >  (] 

where  I0  and  /j  are  modified  Bessel  functions  tabled  in  Abramowitz  and  Stegun  (1964). 
Tables  of  Pfcj^Xj  +  c^'xi  >  c]  are  given  in  Solomon  (1960).  If  tables  of  non-central  chi- 
square  distribution  functions  are  available,  we  may  use  the  exact  expression  that  follows 
where  A  and  B  are  non- centrality  parameters. 

P I5(»xf  +<»  X?  >  I) 

=  P[xl.A  <B1-  P[xl,B  <  A\ 


where 


X  =  ^(3-5i, 
B  =  i(3  +  5>). 


Now  for  r  =  5, 


P  [6x2  +  Xl  >  *] 

=  P  lx?  >  <]  +  p  [x2  <  I]  <x>( \)’p  [*2  >  j 


5*1  /6 


=  P[x2>t]+.l(^)3e^  + 


+  «tt(.6912  +  .144f)P 


fxl  <  t!  • 


The  asymptotic  distribution  of  f  is  summarised  in  the  table  below  for  small  values  of 


r: 


number  of  items 
ced 


Q  =  asymptotic 
distribution  of 


P[Q>t\ 


3 


4 


(nrn  + 


4x!  +  Xi 

S^xS  +  ^xS 


P[5(1)x?  +  «x?  >  I]  +  e~  3t^.{I0{.2t) 
+1.5/1(.2<)} 
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20.  'ABSTRACT 


Exact  expressions  for  the  distribution  function  of  a  random  variable  of  the 

form 

2  2 

c,X„  +  c,x„ 

i  m  /  n 

2  2 

are  given  where  Xm  and  XR  are  independent  chi-square  random  variables  with 
m  and  n  degrees  of  freedom  respectively.  (The  positive  c^  are  distinct.) 

In  particular,  the  exact  asymptotic  function  for  the  average  Kendall  tau  statis¬ 
tic  is  written  as  a  function  of  tables  of  Solomon  (1960)  and  some  found  in 
Abramowitz  and  Stegun's  Handbook  of  Mathematical  Functions. 


